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Starting with a microscopic model based on the Kohn-Luttinger Hamiltonian and kinetic p-d ex- 
change combined with Boltzmann formula for conductivity we identify the scattering from magnetic 
Mn combined with the strong spin-orbit interaction of the GaAs valence band as the dominant mech- 
anism of the anisotropic magnetoresistance (AMR) in (Ga,Mn)As. This fact allows to construct a 
simple analytical model of the AMR consisting of two heavy-hole bands whose charge carriers are 
scattered on the impurity potential of the Mn atoms. The model predicts the correct sign of the 
AMR (resistivity parallel to magnetization is smaller than perpendicular to magnetization) and 
identifies its origin arising from the destructive interference between electric and magnetic part of 
the scattering potential of magnetic ionized Mn acceptors when the carriers move parallel to the 
magnetization. 

PACS numbers: 71.70.Ej,72.25.Rb,75.47.-m 



I. INTRODUCTION 

Although it has been known for 150 years that electric 
resistance of a magnetic metal depends on the direction 
of magnetization^ the origin of such dependence is often 
explained only vaguely as an interplay of spin-orbit inter- 
action (SOI) and magnetization. Conceptual questions 
around this phenomenon, the anisotropic magnetoresis- 
tance (AMR), remain open and relate to the quest for 
its detailed mechanism, its sign and specific magnitude. 
Our ability to control the AMR by material design, with 
potential impact on new electronic devices^ would be 
improved if we had answers more specific but neverthe- 
less still more universal than "black-box-like" modelling 
of AMR in each and every dcemablc system. 

Several factors obstruct a clearer insight into the phe- 
nomenon: there are many electronic bands crossing the 
Fermi level in most materials and the AMR of a crys- 
talline material has various contributions of different 
symmetries. Ab initio calculations performed in FeNi 3 - 
and FeCo4 disordered alloys agree reasonably well with 
experimentally determined AMR but they do not allow 
for any detailed conclusions about its mechanisms. On 
the other hand, the model^ of current-carrying s-states 
scattered to spin-orbit-coupled c?-states provides a rel- 
atively transparent picture of the AMR£ but requires 
the fitting of one or more phenomcnological parameters 
and even then a clear-cut correspondence to ab initio 
results for ferromagnetic transition metals could not be 
established^ 

Diluted magnetic semiconductors, and (Ga,Mn)As in 
particular, offer a promising system in which these is- 
sues become simplified £ Fermi level lies close to the 
top of the valence band so that only few bands are in- 
volved in transport, and in addition, their SOI is strong. 
Moreover, experiments done so far show that the non- 
crystalline component of the AMR^ arising from the 



breaking of the symmetry by choosing a specific cur- 
rent direction, outweighs the crystalline components in 
most of the metallic highly Mn-doped materials. In at- 
tempting to describe the AMR in such system, we can 
begin with a model isotropic but still spin-orbit coupled 
band structure and add the effect of magnetization in the 
three possible distinct ways, as sketched in Fig. Q] Ei- 
ther (a) the magnetization induces a magnetotransport 
anisotropy via the SOI already at the level of group ve- 
locities of the exchange-split Fermi surfaces^ (b) it may 
enter via anisotropic scattering of the unpolarized spin- 
orbit coupled carriers from polarized magnetic impuri- 
ties, or (c) via anisotropic scattering of partially polarized 
carriers which does not require a magnetic character of 
the scatterers. We point out that the mechanisms (a) and 
(c) represent a situation where both fundamental ingrc- 
diences of the AMR (SOI and magnetization) are present 
in the same states of the band structure. The SOI is nec- 
essary for AMR to occur but, at the same time, it weak- 
ens the effect of magnetization so that weaker AMR may 
be expected whenever the mechanisms (a) or (c) dom- 
inate. On the other hand, in mechanism (b), the SOI 




FIG. 1: The AMR can originate from three distinct mech- 
anisms combining magnetization (pointing to the right 
in this sketch) and SOI, that break the isotropy: (a) 
anisotropic Fermi velocities (arrows) along the Fermi surface, 
or anisotropic relaxation rates due to (b) unpolarized bands 
(represented by the indicated isotropic spin texture) scattered 
by anisotropic impurities or (c) partially polarized bands scat- 
tered by isotropic impurities. 
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in an unpolarizcd carrier band can be strong while the 
magnetization of the impurities remains at 100%. Con- 
sequently, very large AMR can arise if this mechanism is 
important^ 

We show in this paper that metallic (Ga,Mn)As is a 
favourable system for the purposes of studying AMR. Not 
only because of its relatively simple (effective) Hamil- 
tonian (described in Sec. [TTJ) and the dominance of the 
AMR mechanism (b), but also because of the way the 
AMR model can be simplified (as shown in Sec. 
down to analytical formulae revealing the basic AMR 
trends (see Sec. IIV[) . This analysis is our main result 
together with the detailed explanation of the AMR sign 
in (Ga,Mn)As (resistance parallel to magnetization is 
smaller than perpendicular to magnetization) which is 
observed experimentall y 8 ' 11 ! 12 ! 13 ! 14 ' 15 ! 16 ! 17 and is oppo- 
site to most magnetic metals . 18 ! 19 The results in Sec. IIVI 
include analytically evaluated anisotropic conductivity 
on several levels of model complexity, and the most 
simplified model allows to clearly identify the physi- 
cal mechanism that determines the sign of the AMR 
in (Ga,Mn)As. Our approach^ is based on the relax- 
ation time approximation (RTA). It is desirable to put 
the present results into more precise terms by exactly 
solving the Boltzmann equation in its full integral form as 
the authors did for the simpler Rashba system recently^ 
However, we have also shown in detail that the sign and 
the qualitative behaviour of AMR in model systems can 
be understood using the RTA 3 The explanation of the 
negative AMR (resistance parallel to magnetization is 
smaller than perpendicular to magnetization) provided 
here is robust against the choice of approximation. 



II. BASIC MODEL OF AMR IN METALLIC 
(Ga,Mn)As 

Three principal ingredients, described in the following 
Subsections (A,B,C), are necessary to model the con- 
ductivity and its magnetic anisotropy: (A) The band 
structure yielding the spectrum and wave functions, (B) 
the scattering mechanism, and (C) a transport formalism 
which combines the former two and produces the conduc- 
tivity tensor. Given that we base our approach to (C) on 
relaxation-time approximate solution to the semiclassical 
Boltzmann equation, we basically need the Fermi veloci- 
ties derived from the band dispersions, and the relaxation 
times calculated from the spectrum, wave functions and 
the relevant form of the impurity potential. 



A. Virtual-crystal kinetic-exchange model of 
(Ga,Mn)As bands 

The valence-band kinetic-exchange model of 
(Ga,Mn)As with metallic conductivities is an estab- 
lished qualitative and often semiquantitative theoretical 
approach £22- The description is based on the canon- 



ical Sctiricffcr- Wolff transformation of the Anderson 
Hamiltonian- 2 ^. which for (Ga,Mn)As replaces hybridiza- 
tion of Mn d orbitals with As and Ga sp orbitals by 
an effective spin-spin interaction of (L = 0; S = 5/2) 
local moments with host valence-band states. This step 
proves essential to effectively separate the different AMR 
mechanisms (a,b,c), symbolized in Fig. [TJ because - 
except for the spin-spin interaction which will be treated 
as we describe below — it completely detaches the Mn 
states from the spin-orbit coupled host-valence-band 
states. These valence-band states are conveniently 
parametrized by the Luttinger parameters 71,72,73 and 
spin-orbit splitting Aso in the six-band Kohn-Luttingcr 
Hamiltonia n 24 ' 25 Hkl- The local interaction between 
Mn magnetic moments Si (located at Rj) and valence 
hole spins s (at r), being at the root of the carrier- 
mediated fcrromagnctism in (Ga,Mn)As, is the kinetic 
exchange and it is described by single parameteri^ 6 - J p d- 
In order to model the band structure of (Ga,Mn)As 
including disorder electrical potential V associated with 
the Mn magnetic moments, we treat the Hamiltonian 

H = H KL + V dls = (1) 

H KL + J P d Sl ■ s S( - r - R ') + Yl V ( r _ R/ ) 
1 1 

by the virtual-crystal mean-field^ 7 - approximation, 
whence we get the single-particle Hamiltonian (in mo- 
mentum representation) of the (Ga,Mn)As valence band: 

H = H KL + he M ■ s . (2) 

Here, &m stands for the unit vector in the direction of 
the mean-field magnetization, h = JpdNMnSM.ni and the 
magnetic moment of Mn is Sjytn = 5/2- In this paper, we 
will only consider substitutional Mn with volume density 
A^Mn as in optimally annealed samples,— and assume zero 
temperature. In the band structure model, we thus dis- 
regard the randomness in the Mn distribution over the 
crystal and the ensuing spatial inhomogeneity of the ex- 
change interaction, and also we completely ignore the 
disorder defined by the electrical potential V in Eq. (TT]) 
of every single substitutional Mn which is an ionized ac- 
ceptor. Within this approximation, the effect of the Mn 
atoms present in the crystal is reduced only to the ef- 
fective Zeeman-likc term in Eq. @ due to the kinetic 
exchange of the valence holes with the Mn d-statcs. Ex- 
plicit form of the k-dependent 6x6 matrix Hkl in a 
convenient basis is given e.g. by Eq. (A8) of Ref. [27l . 

As we are aiming at a simple model of the non- 
crystalline AMR component only, we will treat Hkl 
in the spherical approximation, implemented by setting 
72, 73 to their average valued In this approximation the 
dispersion of all six valence bands becomes isotropic in 
the absence of the kinetic-exchange field. The 6x6 
Hamiltonian t^fy can be diagonalized numerically and 
provide the valence bands E n (k) of (Ga,Mn)As which 
are split by the exchange field h. The index n labels the 
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two heavy hole bands (n = 1,2), two light hole bands 
(n = 3,4), both of the symmetry and total angular 
momentum J = 3/2 in the T-point, and two split-off 
bands (n = 5,6) with the I\ symmetry and J = 1/2 
in the T-point. Note that spin is not a good quantum 
number owing to the presence of SOI. Expectation value 
of spin along any of the Fermi surfaces can be visual- 
ized as a spin texture rather than having separate spin 
up and spin down bands. An example in Fig. [TJb) that 
corresponds^ to the n = 1,2 bands of Hanhltonian @ 
with h — * 0, shows that for each k there are two states 
with opposite spin whose direction, however, depends on 
k, contrary to systems without SOI. 



where we use^°- 9 VV > , the angle subtended by the velocities 
v„(k) and v„/(k') to take into account v which need not 
be parallel to k in case the combined effect of the SOI and 
magnetization distorts the Fermi surfaces as suggested by 
the sketch in Fig. HJa). 

Substitutional Mn act as acceptors and their mag- 
netic moments participate in the ferromagnetic order 
of (Ga,Mn)As. Acknowledging the magnetic and non- 
magnetic part of Vdis, we take 

M% = (z k , n ,\M B + M c \z kn ) (4) 

for the scattering matrix elements between two eigen- 
states of the Hanhltonian (J5J. In the six-band notation 
of Eq. ([2]), the magnetic part of a single Mn impurity 
scattering operator is 

M = JpdSMn&M • s , (5) 

corresponding to the second term in Eq. (P). Explicit 
form of the spin 6x6 matrices s is again given in Ref . I27I 
The non-magnetic part M c describes screened Coulomb 
attraction of the valence holes to the ionized acceptors 
and we therefore take 

M c = y(|k-k'|)i, ^) = -|_L_, (6 ) 

6 q -t q TF 

where 1 denotes a 6 x 6 unity matrix, e is the host 
semiconductor dielectric constant, qpp = \/z 2 g/£ the 
Thomas-Fermi screening wavevector ) 20 i 31 and g the den- 
sity of states at the Fermi level. 

It is important that the two scattering operators (|5l6p 
add up 'coherently' in Eq. ((4]). The 'incoherent' sum 
l^wl 2 = \{zuy\M B \z kn }\ 2 + K^'lM^kOI 2 would 
describe a physically different situation with two distinct 
types of scatterers, magnetic and non-magnetic ones. 



B. Scattering on random Mn impurities 



In order to get hnite conductivity at zero tempera- 
ture, we need to go beyond the virtual-crystal concept of 
Eq. ([2]). We follow Ref. [13 and use the Fermi golden 
rule (or first order Born approximation treatment of 
Vdi S ) a s the simplest model of scattering to calculate the 
transport scattering rates r„.k of the Bloch states from 
Eq. ©: 



(3) 

I 

Such incoherent sum, with appropriately defined scat- 
tering operators, was used earlier— to describe more re- 
alistic (Ga,Mn)As systems that contain interstitial Mn 
atoms or As antisites in addition to the substitutional 
Mn. 

To summarize our model description of substitutional 
Mn impurities in GaAs, the Mn atoms in Gai-^Mn^As 
enter our model at three different places: (i) As acceptors 
and in the absence of other dopants they determine the 
Fermi level Ep and therefore the density of states, (ii) 
On the virtual crystal approximation (VCA) level, they 
cause the ferromagnetic-exchange splitting of the hole 
bands, and (iii) because of the random distribution in 
the lattice, the Mn impurities also cause scattering. The 
essential feature of the Mn impurity potentials for the 
AMR is that they contain components which are propor- 
tional to the Mn local moments and that these moments 
arc ordered in the ferromagnetic state, as expressed in 
Eq. (|5|). We stress that considering (ii) and (iii) simul- 
taneously leads to only a small "double-counting" error 
in the description of the effect of the Mn-related impu- 
rity potential. In terms of the VCA, we assume in (ii) 
that each site on the cation (Ga) sublattice of the host 
semiconductor is occupied by a mixture of x Mn and 
1 — x Ga. This yields an effective mean potential which 
shares the full periodicity of the host zinc-blende lattice. 
Strictly speaking, the scattering potential of randomly 
distributed Mn on the cation sublattice should be de- 
scribed as the difference between the full impurity poten- 
tial due to Mn and the above VCA potential. Similarly 
the remaining sites occupied by Ga should be described 
by the difference between the Ga potential and the VCA 
potential. Ignoring the latter difference and taking the 
full Mn impurity potential for sites occupied by Mn when 
describing scattering in (iii) is therefore not a precise pro- 
cedure but it introduces only a small error for Mn dopings 
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not exceeding several per cent. 



C. Conductivity of (Ga,Mn)As in the relaxation 
time approximation 

We now calculate the conductivity tensor using the 
semiclassical formula based on the Boltzmann transport 
equation^ 

°V = e 'E J j^(^)-^ n (k)vi(k)S(E F E(k)) . 

We assume zero temperature hence the conductivity 
is determined exclusively by states on the Fermi level 
Ep. The Fermi velocities are calculated as v„(k) = 
(l/7i)Vk-EW- The conductivity tensor depends on the di- 
rection of magnetisation 8m through Eqs. ^ty and (J5|), 
that is owing to the combined effect of magnetisation 
and SOI. Generally, the tensor ([7]) may be non-diagonal 
and the resistivity tensor pij is its inverse. 

Here we consider a 'pure-AMR-configuration' where 
<Jij is symmetric, i.e. free of any antisymmetric Hall 
components. It can be envisaged as a Hall-bar device 
fabricated from a thin (Ga,Mn)As film with an arbitrary 
in-plane magnetization. Experimentally, the magnetiza- 
tion will be controlled by a weak magnetic field whose 
direct effect on the AMR will be neglected. The lon- 
gitudinal and transverse voltage drops are proportional 
to p xx and p xy ; spherical approximation (see Sec. Ill Ap 
makes both the orientation of the film and of the Hall 
bar device with respect to the crystallographic axes ir- 
relevant. In an out-of-plane configuration, the resistiv- 
ity acquires an antisymmetric anomalous Hall compo- 
nent, see Ref. for a comprehensive review. For the 
in-plane configuration considered here, the resistivity is 
completely symmetric, p xy (-8 M ) = Pxy{e. M ) = P V x{8 M ), 
and obeys PxyW/Pav = Ci sin 2(f) where <fi is the an- 
gle between magnetization and current (Hall bar device) 
direction. This is a result of a general analysis 33 of 
isotropic systems with symmetry broken by the current 
flow. According to this analysis, the diagonal resistivity, 
Pxx{4>) = Pati(l + Cj cos 20), carries the same informa- 
tion about AMR as p xy (</>), which is concentrated into 
the non-crystalline AMR coefficient Cj [p av is the angu- 
lar average of p xx (4>)}- 

We use the following definition of the AMR: 

AMR = -2 °"" - — = 2 Pil - P± , (8) 
<T|| + crj_ p|| + pj_ 

where o\\ and aj_ (p\\ and p±) are the longitudinal con- 
ductivities (resistivities) for current parallel and perpen- 
dicular to the magnetization, respectively. Note that the 
off-diagonal resistivities vanish for <j> = or 90°. For 
practical purposes, we can set ej« in the x-direction and 
£T|| = a xx , <j± = <jyy as long as we stay with the spherical 
approximation. The AMR of Eq. © equals to 2C/. 



III. NON-CRYSTALLINE AMR OF HEAVY 
HOLES IN THE SPHERICAL APPROXIMATION 

In this Section, we first show that within our descrip- 
tion of metallic (Ga,Mn)As samples, the AMR trends are 
governed by the mechanism sketched in Fig.QJb). Then 
we proceed to showing that only heavy hole bands need 
to be considered when analyzing the basic AMR charac- 
teristics in (Ga,Mn)As on a qualitative level. 



A. Origins of anisotropy 

Let us consider how the conductivity in Eq. |(7J) can be- 
come magnetization-dependent. In agreement with the 
intuitive analysis of Fig. [TJ magnetization direction 8m 
can enter Eq. ([7]) cither via (a) the group velocity com- 
ponents v l n or (b,c) the scattering rates r„k. Considering 
Eqs. (|3|4| . the scattering rates may depend on 8 m either 
through (b) the scattering operator (Mb + Mc in our 
specific case) or (c) the wavefunctions |zkn) and/or ener- 
gies E n (k) and density of states of the carrier bands. The 
last mechanism, for example, lies at the heart of the s- 
d model of AMR in transition metals^ where isotropic 
and spin-independent scattering operators have been 
assumedi 18 ' 19 The anisotropy arises due to the compe- 
tition of SOI and magnetization which splits the five d 
states according to Im, their angular momentum projec- 
tion along 8m- These states play the role of j^k'n') in 
Eq. (|4|), and because of their ^/-dependent spatial form 
they cause ejv/-dcpcndent scattering rates r„k in the s- 
states that carry the current. Plugged back into Eq. (J3)), 
these anisotropic scattering rates may lead to different 
a xx for 8m parallel and perpendicular to the a;-direction, 
i.e. £7|| ^ ax- 

In (Ga,Mn)As, we are going to take advantage of the 
tunability that the effective model outlined in Sec. |TT] of- 
fers: we will switch the particular mechanisms (a,b,c) on 
and off to see how important they are for the total AMR. 
We use the band-structure model as described in Sec. Ill Al 
(termed 'full spherical') as a reference. Within this full 
spherical model, the calculated AMR as a function of Mn 
doping x is negative in the considered range between 2 
and 10% and its magnitude reaches a clear maximum as 
shown by the middle curve in the left panel of Fig. [2] (the 
maximum is related to the competition between the elec- 
tric and magnetic parts of the scattering operator as we 
explain in Sec. IIV[) . To see the effect of the mechanism 
(b) alone, we set h — in Eq. ([2]) but leave the scatterer 
anisotropy unchanged by keeping nonzero J p d in Eq. (|5j). 
The bottom curve in the left panel of Fig.[2]demonstrates 
that the AMR quantitatively changes within a factor of 
two but its overall form remains the same. 

On the other hand, the result alters dramatically when 
we switch off the anisotropy in the scattering operator 
[mechanism (b)] or the anisotropy in relaxation rates as 
a whole [mechanisms (b) and (c) together]. The former 
is accomplished by setting Mb = in Eq. the latter 




FIG. 2: Left: study of the influence of the three mechanisms (a,b,c) sketched in Fig.[T]on the total AMR. The model in spherical 
approximation as of Sec. |IT1 serves as a reference ('full spheric'). Note that the upper two curves are upscaled by a factor of 
ten. Right: subsequent approximations (see text) aiming towards an analytically solvable model. Note that (B) and (C) is the 
same as 'full spheric' and (b) on the left panel, respectively, while (D) and (E) correspond to gradual decoupling of the light 
holes. 



is done by replacing T n ^ by a constant whose value is ir- 
relevant because it cancels out in Eq. ([8]). In both cases, 
we obtain AMR that is more than an order of magnitude 
smaller than for the full spherical model, see Fig. [2j This 
result indicates that mechanism (b) is a crucial part of 
the AMR model of metallic (Ga,Mn)As, the other mech- 
anisms (a,c) provide only quantitative corrections and 
when kept alone without mechanism (b), they produce 
negligible AMR. 

B. Heavy holes 

The next simplification of the model we make in or- 
der to provide a simple physical picture of the AMR in 
(Ga,Mn)As, is to neglect the light holes. We can accom- 
plish this in two steps: we first discard the current carried 
by the light holes, i.e. sum in Eq. ([7]) over n = 1, 2 only, 
and then also disable the scattering from the heavy hole 
to light hole bands, i.e. sum in Eq. ([3]) over n' = 1,2 
only. Numerical calculation again shows that this pro- 
cedure does not alter the qualitative behaviour of the 
AMR. All levels of approximations are summarized in 
the right panel of Fig. [2j For completeness, we start with 
the 'full spheric' reference (curve labelled by 'B') as on 
the left panel and proceed to suppressing the band polar- 
ization by putting h = in Eq. ^ which yields the curve 
'C. Omission of light-holes from current-carrying states 
in the transport equation but not from the final states in 



the scattering matrix elements produces the data labelled 
by 'D', and the completely heavy- hole-only model (six- 
band model where the other four bands are disregarded) 
is denoted by 'E\ 

For further studies of this model it may be interesting 
that the difference between curves 'C and 'D' is remark- 
ably small, in other words, the anisotropy of the light hole 
transport is almost identical to that of the heavy holes, 
provided we have set ft = in Eq. (J5J. However, the 
main conclusion of this Section is that, within the stud- 
ied range of Mn doping, the AMR is determined by the 
anisotropy of the relaxation rates of heavy holes induced 
by the Mn electro-magnetic scatterers. This allows us to 
derive an approximate analytical formula for the AMR 
in (Ga,Mn)As which we discuss in the following section. 



IV. QUALITATIVE ANALYTICAL RESULTS 
FOR AMR IN (Ga,Mn)As 

We first provide analytical expressions for the AMR 
corresponding to the curve 'E' in Fig.[2]with an additional 
approximation that the current I is proportional only to 
the transport life-time of carriers with v„ || I. 

Using the explicit form of the heavy hole wavefunctions 
given in Ref. [29|, the transport scattering rates of Eq. ([3]) 
corresponding to M B +M C of Eqs. (|5l6p can be evaluated 
as 
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r±(o) 



6 + 1 



T ( 12 + 186 + — - - (1 + 36) 2 In — - J T a [ 4 + 66 + 66 2 - (6 + 1)(1 + 36 2 ) In 
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2A: 2 , 



(9) 



where r±(0) are the scattering rates of the two heavy-hole bands with <f) equal to the angle between k and &m- Next, 
we get 



186 



(l + 36) 2 ln- 



1 



(10) 



We used shorthands a = — (e 2 /2ek F )/J p dSMn with &f 
denoting the common Fermi wavevector of the heavy 
holes, and 6 = —(1 + q FF /k 2 F ). 
Conductivities are evaluated as 
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-(0) 
1 



+ 



r_(o) 



r + (H r_(i7r) r+fa)' 

and yield AMR which is up to an overall factor of w 2 
the same as if we performed the complete k-integration in 
Eq. ([7]), see Fig. [3Ja). The approximation of taking into 
account only states with k || I is therefore qualitatively 
valid. 



Short range scatterers 



Results of Eqs. (|9|10|) are analytical but the formula 
for AMR is rather complicated. We can attain a clear 
insight into how the observed AMR trends arise if we 
further simplify the model. We replace M B + M c of 
Eqs. (|5I6[) by another scattering operator 



M s - r - tx e M -s + al 



(13) 



That is, we assume k-independent electric part of the 
scattering potential. 

The transport scattering rates of Eq. ((3J again depend 
only on the angle <f>, 

r±(0) oc i cos 2 4> ± acos0 + a 2 + — , (14) 

and specifically for states with k parallel and perpendic- 
ular to &m, we get 

r ± (0)(x(a±i) 2 , r ± (i7r)cxa 2 + ^. (15) 



Conductivities evaluated using Eqs. (|ll|12p simplify to 

1 1 



(7|| CX 



er_L oc 



(a + 1)> T (a 



r+ (W) a2 + i' 



(16) 
(17) 



p[nm ] 
0.44 0.88 1.32 1.76 



p[nm ] 
7.4 1.76 0.44 



0.13 




FIG. 3: (a) The analytical model of AMR based on 
Eqs. (|9I10|I . Curves 'B', 'E' are taken from the right panel 
of Fig. [2] Estimates of the hole concentrations shown on the 
upper axis assume one hole per Mn and neglect the light- 
hole bands, (b) The AMR as a function of the effective 
Coulomb scattering strength given by Eq. (|19[l . The 'limited- 
integration' model (only states with k parallel to the current) 
is represented by Eq. QS) , the integrated model where all 
states contribute to the current is described after Eq. (I20I) . 



which gives, using Eq. (J5]), our previous result^ 2 ^ 



AMR 



20a 2 



24a 4 - 2a 2 



1 



(18) 



In order to link this result, plotted in Fig. [3fb), to the 
previous one that is based on the full form of the Coulomb 
scattering operator given by Eq. ([6]), we need to focus 
on the parameter a. It represents the effective strength 
of the electric part relative to the magnetic part of the 
scattering potential of the Mn ions. We can estimate a 
as an average over the Fermi surface of the more realistic 
V(|k-k'|) from Eq. ©, 



(V) FS _ e 2 /e 

JpdSMn JpdSMn 



1 

4fcf 



In 



QtfJ 



(19) 



Explicitly, (V) FS = {^k 2 F )- 1 j FS d 2 k'V(\k~ k'|) with k 
fixed to an arbitrary Fermi wave vector, |k| = kp. The 
integral is taken over the Fermi surface. 

We thus find that the relative strength of the Coulomb 
scattering increases with decreasing hole density, p w 
kp/(3ir 2 ), due to the screened long-range nature of the 
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Coulomb interaction which contrasts the short range of 
the magnetic scattering. Namely, the momentum trans- 
fers q scale down with kp and we increasingly get to 
feel the singularity of V(q) at q = even though it is 
rounded by the screening, see Eq. ([6]) . Assuming a value 
p = 0.5 nm -3 which is realistic in (Ga,Mn)As, we have 
kp ~ 2.5 nm _1 , qpp * 1-6 nm _1 , and a m 1.0. 

Equation (|19[) provides the link between a and p which 
is determined by x, i.e. between horizontal axes of 
Figs. [3£a) and (b). Qualitatively, increasing x and thus 
also increasing kp implies decreasing a because the scat- 
tering rate due to the magnetic part scales oc kp while 
the scattering rate due to the Coulomb scattering given 
by M c of Eq. ([6|) grows considerably slower. If we es- 
timate V(|k — k'|) by V(2kp), the growth will be only 
oc k F I (kpa 2 + kpaj ir) where a is the effective Bohr radius 
in GaAs. Behaviour of the more sophisticated estimate 
used in Eq. (|19[) is qualitatively the same. 

The simple form of scattering rates in Eq. (TT4")) al- 
lows to analytically perform the full integration over k 
in Eq. ([7]) and to obtain more precise results for conduc- 
tivity. For ctii = <t^ + o~7 we get 



FS 



cPk cos 2 i 



(20) 



The proportionality factor is e 2 times squared Fermi ve- 
locity times the density of states at the Fermi level. Con- 
ductivities of the two bands are equal, so that 



12 - 36a In 



(21) 



24a 2 



: arcsinh 



1 



^|6a 2 - 1|/18 
l« 2 -il 



a > 7 
6 



and the same result with arcsinh replaced by arcsin ap- 
plies for a 2 < i. 

In order to evaluate <tj_, we employ a straightforward 
identity + 2crJ= = T± with 



2V18 



i arcsinh 



V|6« 2 - 1|/18 



o- 



a > - 
6 



and again with arcsin for a 2 < i, which can be derived 
by inserting cos 2 0=1- sin 2 <p into Eq. ([20]). The AMR 
evaluated using Eq. JS]) is shown as the solid curve in 
Fig. lib). 

Both models with simplified scattering operator M s r - 
exhibit qualitatively the same behaviour of AMR as a 
function of a: positive AMR for a close to zero, sign 
change to negative AMR already for a small value of a, 
maximum AMR magnitude at a = 1/2, and vanishing 
AMR for ol — ► oo. A comparison between Figs.OHa) and 
(b) reveals that the outstanding feature of the AMR at 
| a | = 1/2 makes its way2& up to the full spherical model 
where it is becomes broadened to the wide maximum 
(around x w 4% in Fig. [5]). Based on the model that 



employs M s - r - of Eq. (|13p . we now analyze the origin of 
the very large negative AMR at |a| = 1/2. 

This maximum AMR value follows from the diverg- 
ing conductivity erii given by Eqs. (fTB|l or (|2"T]) which is 
caused by vanishing scattering rate r =F (0) for a = ±1/2. 
The scattering matrix element in Eq. Q vanishes for 
any n' , k', and so does the integral (J7J), when |zk n ) is an 
cigenstate to M s ^ r ^ with eigenvalue zero; for illustrative 
purposes, consider M s ' r ' = s x + a, i.e. % pointing to 
the right in terms of Fig. Q] The absence of scattering 
occurs for a = 1 /2 when k is parallel to the x-direction, 
as evident^ 9 - from the spin texture in Fig. \V[b). Such 
infinite conductivity will of course not occur in realis- 
tic systems; as soon as the dependences on momentum 
transfer of the electric and magnetic parts of the scatter- 
ing potential will not be exactly the same, the eigenstate 
property of \zk n ) is lost. This is the case when we replace 
M s - r - by the original M B + M c of Eqs. (|5l6"f . However, 
as Fig. [3f a) shows, there still remains a maximum in an 
which translates into a maximum of |AMR| as a finger- 
print of the original cry singularity. 

Positive AMR for purely magnetic scattering (a = 0) 
that follows from Eq. (fig)) can also be understood us- 
ing the basic properties of spin-3/2 states.— The idea is 
that |zkn) is an eigenstate (with non-zero eigenvalue) to 
■ s for k oriented along &m thereby allowing scatter- 
ing to a state with — k that contributes strongly to the 
transport relaxation rate of Eq. ((3|) . On the other hand, 
matrix elements of eu ■ s vanis h 10 ' 29 between states with 
k, k' _L 8m and scatterring and resistivity are therefore 
suppressed for the magnetization perpendicular to the 
current. Relation a± > an (AMR > 0) for a = 0, ob- 
tained from the above analysis of the leading scattering 
channels which contribute to transport, is confirmed by 
calculations that even take into account all final states 
k' for the scattering and include the integration over all 
k-states in Eq. 0. The sign of the AMR changes quickly 
from positive to negative when adding the electric com- 
ponent of the scattering potential. In the model of short 
range potentials, summarized by Eq. (|18|) . it occurs at 
a = 1/V20. Recall that for realistic impurities with 
screened long-range Coulomb potential the transition can 
be triggered by changing hole density, as explained in the 
discussion of Fig. O 

At the end of this Section, we recall that all consid- 
erations in this work are based on the RTA. It has been 
shown in Ref . [2l| that the RTA does not provide the exact 
solution to the Boltzmann equation. Even in the simple 
case of a spherical model in (Ga,Mn)As we should solve 
an integral equation rather than merely evaluate inte- 
grals as in Eq. ([3]). Nevertheless, it has also been argued 
in Ref. [l^ that the basic qualitative trends of the AMR 
can be found already on the level of the RTA. In partic- 
ular this applies to the maximum of |AMR| at |a| = 1/2 
whose robustness was verified by comparing the RTA and 
the exact solution of the Boltzmann equation in its inte- 
gral form on model Rashba-Drcssclhaus SOI systems. 
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V. SUMMARY 

Non-crystalline anisotropic magnetoresistance (AMR) 
is governed only by the angle between magnetization and 
current rather than by their orientation with respect to 
crystallographic axes, as opposed to the crystalline AMR. 
The physical origin of the AMR is the combination of 
spin-orbit interaction (SOI) and of the broken symmetry 
due to the presence of magnetization. More specifically, 
three distinct mechanisms may lead to the non-crystalline 
AMR: Fermi surfaces distorted from the spherical shape 
that imply anisotropic Fermi velocities, and anisotropic 
scattering rates, either due to anisotropic wavefunctions 
or due to the anisotropic scatterers. We note that the 
s-d model which is sometimes invoked to qualitatively 
explain the AMR in ferromagnetic transition metals, is 
a variation of the anisotropic wavefunction mechanism. 
Because of the competing effect of the SOI and magneti- 
zation (which are both present only in the low-mobility 
d states) one may expect weak AMR. 

On the other hand, we have shown that (Ga,Mn)As 
constitutes a textbook example of a system where strong 
AMR can be expected since the two agents are mostly 
separated: polarized Mn ions act as anisotropic scatter- 



ers while the current-carrying valence-band states bring 
in the SOI. Quantitatively, the latter are also partly po- 
larized but we have demonstrated that for the typical ex- 
perimental range of Mn dopings in metallic (Ga,Mn)As 
samples [x = 2 ~ 10%), the anisotropic scatterer mech- 
anism is dominant. A simple model which neglects the 
other two mechanisms provides analytical results which 
predict the correct sign of the AMR (resistivity parallel 
to magnetization is smaller than perpendicular to mag- 
netization) and identify its origin — destructive interfer- 
ence between electric and magnetic part of the scattering 
potential (of ionized Mn acceptors) for carriers moving 
parallel to magnetization. 
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NAMASTE (FP7 grant No. 214499) and SemiSpinNet 
projects (FP7 g rant No. 215 368), SWAN-NRI, ONR un- 
der Grant No. Ionr-n0001406l 0122. NSF under Grant 
No. DMR-0547875, and also Praemium Academias. 
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illustrated on Fig. [T] Consider the spin texture in Fig. []Jb) 
as a typical example of SOI effect. Exchange splitting due 



to magnetization oriented horizontally in Fig. [TJ effectively 
acts as a Zeeman energy due to fictitious 'magnetic field' 
B with the same orientation, see Eq. The states with 
wavevector k||B will be shifted upwards (downwards) in 
energy because their spin is parallel (antiparallel) to B 
while the states k 1 B remain intact. This k-anisotropic 
shift in energy can be translated into the splitting of Fermi 
wavevectors of the two originally degenerate bands. 
Note that once we have set h = in Eq. ([2]), the theoretical 
AMR depends on the Mn doping x only through the carrier 
density p[nm ] = 0.22 x x[%], because A^Mn appearing in 
Eq. drops out in Eq. ||8J). This fact should not be taken 
ad absurdum (Nmu — ► 0) because we still assume that the 
substitutional Mn impurities provide the dominant source 
of scattering. Presence of other concurrent types of impu- 
rities (which is important when we compare the calculated 
AMR to experiments^) will make the AMR depend on the 
ratio between their concentration and Nmh- Nevertheless, 
this interplay will not simply obey Matthiessen's rule since 
we deal with anisotropic systems.— 



